import numpy as np

a = np.array([[1., 2.], [3., 4.]])
y = np.array([[5.], [7.]])

print(a)
print(a.T)

dia_num = [1.0, 1.0]
b = np.diag(dia_num)
print(b)
print(b.trace())
print(np.linalg.solve(a, y))
